o 

O 
(N 



Frequency metrology on the 4s ^Si/2 — 4p ^Pi/2 transition in 

data 



40 



Ca+ for a comparison with quasar 



A.L. Wolf,i-2.Bs.A. van den Berg,^ C. Gohle,' E.J. Salumbides,' W. Ubachs,' and K.S.E. Eikema' 

Laser Centre Vrije Universiteit, De Boelelaan 1081, 1081 HV Amsterdam, The Netherlands 
'NMi van Swinden Laboratorium, Thijsseweg 11, 2629 JA Delft, The Netherlands 

(Dated: May 14, 2008) 

High accuracy frequency metrology on the 4s S1/2 — 4p P1/2 transition in calcium ions is performed using 
laser cooled and crystallized ions in a linear Paul trap. Calibration is performed with a frequency comb laser, 
resulting in a transition frequency of / = 755 222 766.2(1.7) MHz. The accuracy presents an improvement of 
more than one order of magnitude, and will facilitate a comparison with quasar data in a search for a possible 
change of the fine structure constant on a cosmological time scale. 



I. INTRODUCTION 
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The detection of a possible variation in the fine structure 
constant a by Webb et at. iM was founded on the so-called 
many-multiplet method. In this study accurate wavelengths 
of absorption lines of a number of atoms and ions, observed 
at high redshift in the line-of-sight of Quasi Stellar Objects 
(QSO or 'quasars'), were compared with the laboratory (zero 
redshift) wavelengths of those atomic resonances. Such an 
analysis may yield a value for Aa/a over time spans of many 
billion years. Improved methods and inclusion of over a hun- 
dred of high-redshift absorption systems now has delivered a 
positive effect with a 5<j significance |2]. The cuiTent oppor- 
tunities for highly accurate wavelength calibrations of astro- 
physical data, provided by echelle-grating spectrometers at- 
tached to very large telescopes, have brought about the need 
for renewed and updated laboratory precision measurements. 
For several spectroscopic lines the absolute wavelength accu- 
racy does not match the level of 10^^, which is feasible in 
astrophysical observations. One of the lines that need a more 
accurate frequency calibration is the 4s ^Si/2 — 4p ^Pi/2 tran- 
sition in Ca+, which is the subject of this paper p']. 

In general, trapped calcium ions are studied both for use 
in very accurate ion clocks 141 1511 and for quantum compu- 
tation Ii^lZJ. Nevertheless, the 4s^Si/2 — 4p^Pi/2 transi- 



tion that is used for Doppler cooling these ions, has been 
given far less attention. For this transition, often the value 
1/ = 25 191.5182 cm^^ is cited [SJ. This frequency is based 
on a private communication with U. Litzen, who recently re- 
calibrated the coiTesponding data yielding an updated value 
of ly ^ 25191.519(0.002) cm^^ |9]. The next best wave- 
length calibration dates back as far as 1956, with a value 
I' = 25 191.51(0.13) cm~^, obtained by Edlen and Risberg 
using classical spectroscopy lUOll . 

In the present study, a linear Paul trap is used for con- 
finement of a cloud of calcium ions, while laser cooling is 
employed to obtain a linewidth equal to twice the natural 
linewidth of the transition. The cold ions are probed using 
a second laser beam of lower intensity. Simultaneously, a part 
of this laser is heterodyned to a frequency comb for absolute 
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II. EXPERIMENTAL SETUP AND PROCEDURES 

A. Ion production and trapping 

The calcium ions are trapped in a linear Paul trap (Fig. [1]) 
consisting of 4 linear electrodes for radial confinement, and 
2 ring electrodes for restriction in the axial direction. These 
elements are suspended and separated by insulating spacers. 
The electrode dimensions enable fluorescence collection over 
a large solid angle of « 0.27r sr. The four linear electrodes 
are placed at a distance tq = 4 mm from the trap axis, and 
the ring electrodes are 17 mm apart. A radio frequency (RF) 
voltage Vb at a frequency il ^ 2tt x 3.3 MHz is applied to the 
electrodes 1 and 4, while the electrodes 2 and 3 are grounded. 
The characteristic stability parameter of the Mathieu equation, 
which describes the motion of an ion in an harmonic RF po- 
tential, is q = 2QVo/mil^rQ |11]. Here Q and m are the 
charge and mass of the ions respectively. A small stability 
para meter (q << 1) is advantageous to reduce RF heating 
II12I [l3ll . However, trap imperfections and micromotion in 
an extended cloud can still lead to substantial heating of the 
confined ions Illll[l4ll . The trap is generally operated with a 
Vb ~ 60 V, corresponding to q = 0.04. The ring electrodes 
that are held at a constant voltage of typically Vdc = 10 V 
create a harmonic potential for axial confinement. 

The trap is installed in a vacuum chamber, which is evacu- 
ated to a pressure of 2 x 10^^ mbar. Neutral calcium atoms 
are evaporated from an oven into the linear Paul trap. In the 
trapping region calcium atoms are resonantly excited to the 
4p ^Pi level by a frequency-doubled CW (continuous wave) 
Ti: sapphire ring laser (Coherent 899) producing 422 nm, and 
subsequently ionized by a frequency tripled Nd: YAG (yttrium 
aluminum garnet) laser operating at 355 nm (Fig. |2]i. The RF 
source for the trap (Agilent 33120A) is resonantly enhanced 
and filtered by a helical resonator (with a loaded Q«140). For 
diagnostic purposes a Channel Electron Multiplier (CEM) and 
accelerating grid were installed underneath the Paul trap. By 
pulsed extraction of the ions from the trap with a variable de- 
lay and detection on the CEM, the ionization rate and trapping 
conditions can be probed and optimized. 




(a) Side view 



(b) Axial view 
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FIG. 1: (Color online) Electrode configuration for the linear trap. 
The trap consists of four linear electrodes (labeled 1-4) for radial 
confinement: A Radio Frequency voltage is put on two diagonally 
facing rods (1,4), while the other two rods are grounded (2,3). A DC 
voltage is applied to the ring electrodes for axial confinement. 
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FIG. 2: (Color online) Schematic design of the laser system 
and linear Paul trap. DL=Diode laser, PMT=Photo-Multiplier 
Tube, 2cij=frequency-doubling setup, 3(j;=frequency-tripling setup, 
PBC=Polarizing Beamsplitter Cube, CEM=Channel Electron Mul- 
tiplier. 



B. Doppler cooling to ion cloud crystallization 



FIG. 3: Partial energy level diagram of Ca 



cused into the trapping region. This further reduces the back- 
ground signal due to scattered light from the trap electrodes, 
which mainly originates from the 397 nm cooling laser. In 
Fig. |4|two excitation spectra of trapped and laser cooled ions 
are shown, obtained by scanning the cooling laser. Depending 
on the cooling efficiency (controlled by the repumper detun- 
ing), the recorded fluorescence follows either a red-shifted, 
distorted Gaussian profile (red curve), or shows a dip fol- 
lowed by a much narrower resonance (blue curve) which is 
typical for crystallization [15]. Such a crystallized ion cloud 
is used for the spectroscopy to minimize Doppler broadening. 
Laser cooling is performed with the cooling laser set at detun- 
ing A/ « 10 MHz from resonance (locked to a wave meter 
(Atos model LM-007) within « 6 MHz), which is approxi- 
mately half the natural linewidth. The achieved temperature 
is typically T « 1 K. This is higher than the Doppler limit of 
T = 0.5 mK, probably due to anharmonicity of the applied 
electric field or other trap imperfections Illll[l4ll . The large 
number of trapped ions (^ 100) magnifies this effect, since 
the Coulomb repulsion forces the ions away from the more 
harmonic trap center The lifetime of the ions in the trap is 
measured to be r « 8 min by monitoring the fluorescence 
from the cooling laser. 



In order to counteract RF heating and reduce the Doppler 
width of the measured transition, the ions are laser cooled on 
the 4s ^Si/2 — 4p^Pi/2 transition (the relevant energy levels 
are shown in Fig. O. The cooling laser is a 397 nm grating 
stabilized diode laser {P — i mW, Toptica DLIOO). Since 
the ions in the excited Ap ^Pi/2 state have a 7% probability of 
decaying to the long-lived 'id^T)^/2 state, an additional diode 
laser at 866 nm (P = 1 mW, Toptica DLIOO) is used for re- 
pumping. 

Excitation of the transition is observed through the emitted 
fluorescence at 397 nm. Using a 25 mm diameter lens with 
a focal length of 25 mm inside the vacuum chamber, 5 % of 
this light is collected. The ion cloud is imaged onto a pin- 
hole to remove background radiation and scattered light from 
the electrodes, and its fluorescence is detected on a photo- 
multiplier tube (PMT, PhiUps XP2020Q). The cooling laser is 
spatially filtered by a single-mode optical fiber (reducing the 
effective cooling laser power to _P = 1 mW), before being fo- 
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FIG. 4: (Color online) Observed fluorescence vs. frequency detuning 
A/ of the cooling laser relative to the 4s ^ S1/2 — 4p ^Pi/2 transition 
in Ca+. Red thin trace: Tions ~ 40 K, blue thick trace: Tions ~ 4 K. 
The sudden decrease in fluorescence in the blue curve is the typical 
signature for crystallization of the ion cloud. 



C. Spectroscopy and absolute frequency calibration 

Because the 4s ^Si/2 — 4p ^Pi/2 transition is both the cool- 
ing and measured transition, the cooling laser is periodically 
switched off using a chopper with a duty cycle of 0.5 at a fre- 
quency of 0.3 kHz. The fluorescence due to the spectroscopy 
laser is recorded during these periods. A diode laser at 794 
nm (Toptica DLIOO) is frequency doubled to the 397 nm to 
probe the transition. To minimize heating of the cloud the 
probe laser intensity is a thousand times weaker than the cool- 
ing laser. After spatial filtering a power of only P « 1 /j,W 
remains, which is focused to 0.5 mm diameter in the trap. The 
induced fluorescence is detected with the imaging system and 
PMT described previously (Fig. |2]). The fluorescence of both 
the cooling laser and spectroscopy are separately recorded, 
using a gating technique synchronized to the chopper in the 
cooling laser beam. The cooling and probe laser are guided 
through the same single mode fiber (Fig. |2l), ensuring good 
alignment of the weak probe beam onto the cold ion cloud, 
while also spatially filtering the probe beam. In addition, by 
monitoring the cooling laser induced fluorescence, ion loss 
from the probed region can be observed during the scan. 

The probe laser is calibrated on an absolute frequency scale 
by referencing it to a frequency comb ifmUTIl . A frequency 
comb laser is characterized by two frequencies: The repetition 
frequency f^ep which determines the mode spacing, and the 
carrier-envelope offset frequency fceo, that gives the comb its 
offset from zero frequency. Both frequencies are locked to a 
Rubidium atomic clock (Stanford Research Systems PRS-10), 
which is referenced to the Global Positioning System (GPS), 
fixing the frequency of then*'* mode to /„ = ztfceo+nxfrep- 
The Allan deviation of the frequency comb modes is deter- 
mined from the stability of the repetition rate relative to the 
clock combined with the stability of the Rubidium frequency 
reference. For the typical averaging time of 1 second, the 
Allan deviation is cr = 0.1 MHz for optical modes at 377 
THz, leading to an uncertainty of 0.2 MHz in the calibration 
at 755 THz. The comb laser operates at a repetition rate of 207 
MHz and with a spectrum of 54 nm FWHM around 825 nm. 

An interference beat note is generated between the spec- 
troscopy laser and the comb laser modes, by overlapping the 
near-infrared fundamental output of the spectroscopy laser 
with the output of the frequency comb laser. The comb light 
is spectrally filtered using a grating and curved mirror in a 2f ' 
configuration (Fig. |5]l to reduce the amount of comb modes 
and improve the signal to noise ratio of the recorded beat note 
with the spectroscopy laser. An avalanche photodiode (Menlo 
Systems APD210) detects the 40 dB (bandwidth 400 kHz) 
beat note /;,„ between the two lasers, which is monitored us- 
ing a frequency counter (Agilent 53132A). 

The laser frequency fi follows from the locked frequencies 
fceo, frep, and the measured beat frequency /;,„ through the 
equation: 
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FIG. 5: Schematic view of the setup for heterodyning the CW diode 
laser with the frequency comb. The CW laser is overlapped with the 
comb on a beamsplitter (BC). A grating-based spectral filter is used 
to reduce the number of comb modes on the detector. For clarity a 
lens is drawn, in practice a curved miffor is used. APD=Avalanche 
Photodiode, ph=pinhole, A/2=haIf-wave plate. 



was repeated at a different locking frequency of the caiTier en- 
velope offset frequency in order to determine the sign of fceo- 
In order to calibrate the laser frequency the mode number n 
still has to be determined. The most accurate measurement 
up to date for the 4s^Si/2 — 4p^Pi/2 transition in Ca+ has 
an uncertainty of 60 MHz, which complicates assigning the 
mode number n with confidence. This is solved by repeating 
the measurement at two different comb repetition frequencies, 
sufficiently far apart to distinguish between the adjacent mode 
numbers. As the transition frequency is independent of the 
used comb repetition rate, the measurements coincide for a 
particular mode number n, from which the absolute frequency 
can be determined (see below). 
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The sign of /b„ can be simply deduced from the known 
scanning direction of the laser. The scan over the transition 



FIG. 6: (Color online) An example of the measured fluorescence 
spectrum for the 4s ^Si/2 — 4p ^Pi/2 transition in Ca"*" (black dots) 
and the corresponding Voigt fit (red line), corrected for the trap loss 
during the scan. 



D. Results 

A typical scan of the spectroscopy laser takes about 2 min- 
utes, which is substantial compared to the trap lifetime of « 8 
min. In order to minimize distortion of the line profile, the 
spectroscopy signal is corrected for this ion loss. Since the 
cooling is insufficient to reach the natural linewidth of the 
transition, a Voigt profile, i.e. a convolution of a Gaussian and 
a Lorentzian is fitted to the data. The width of the Lorentzian 
is fixed to 22.4 MHz 1 18], the natural linewidth of the tran- 
sition. The width of the Gaussian component is fitted, and 
typically on the order of 20 MHz. From this fitted value an 
ion cloud temperature of T w 1 K is deduced (Fig. |6j. 

The mode number is determined by measuring the transi- 
tion frequency using different frequency comb repetition rate 



Of/r 



189 MHz and /,. 



207 MHz. This results in two 



sets of possible transition frequencies (Fig. [Tj. The measure- 
ments clearly coincide at a value of / = 755 222 766 MHz, 
which is near the literature value of / = 755 222 740(60) 
MHz. 

Several effects that can introduce systematic and statisti- 
cal errors on the transition frequency have been investigated. 
These include the Stark shifts caused by both the light field 
and the Paul trap, the stability of the comb modes and the Zee- 
man shift caused by the earth magnetic field. The respective 
shifts and statistical uncertainties are tabulated in Table|Il 
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FIG. 7: (Color online) Absolute calibration of the 4s S1/2 — 
4p^Pi/2 transition in Ca'^ is performed by finding the coincidence 
of measurements with two different frequency comb repetion rates 
of 189 MHz (circles) and 207 MHz (triangles). The depicted tran- 
sition frequency / = 2 x /; , since the laser frequency is measured 
in the infrared, after which the laser is frequency doubled for the 
spectroscopy. Consequently the spacing between adjacent points in 
the picture is 2 x frep- The upper graph shows the frequency de- 
pendence on the choice of mode number. The lower graph shows 
the frequency difference between the points in the upper graph, with 
A/ = / - 755 222 766 MHz-2 x n' x 198 MHz, where n' is the 
relative mode number (error bars are smaller than the point size). 
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FIG. 8: (Color online) Transition frequency for the 4s S1/2 — 
4p ^Pi/2 transition in Ca+ with respect to an applied magnetic field 
Bappiied- Three different polarizations were used (see text): right- 
circular (cr^, green triangles), left-circular (a' , red circles) and lin- 
early polarized light {a^ + cr, black squares). 



In order to probe the Zeeman shift, a magnetic field is gen- 
erated in the direction of the laser beam using Helmholtz coils 
(9 windings, diameter 2r — 19 cm). For magnetic fields in 
this direction, a deviation from linear spectroscopy laser po- 
larization can lead to an unequal excitation of the Am = ±1 
transitions. The Zeeman shift in the involved energy lev- 
els {AE = gjmjiJLBB) results in a magnetic field depen- 
dent shift in the transition frequency of ±19kHz//iT for the 
Am = ±1 transitions. For perpendicular magnetic field com- 
ponents no shifts of the line center, but only broadening can 
be expected to first order. These are therefore not considered. 

The line center of the transition as a function of the ap- 
plied external fields was recorded for three different polariza- 
tions: right- and left-circular ((t+ and a^ respectively) and 
linear {— cr+ + a~). The result is shown in Fig. |8] Lin- 
ear fits to the data show Zeeman coefficients of 19(l)kHz//xT 
and 15(l)kHz//^T for the cr+ and a^ components. The devi- 
ation of the latter from the calculated value can be attributed 
to an inaccurate control of the polarization in the interaction 
region, due to depolarization by the Brewster angled vacuum 
windows. For excitation with linear polarization, a first order 
Zeeman coefficient of 0.0 kHz//iT is expected, and it was cor- 
respondingly measured to be 0.2(1.0) kHz//^T. Therefore the 
final calibration of the 4s ^Si/2 — 4p ^Pi/2 transition is based 
on an average over the measurements with linear polarization. 
In this case, for all applied magnetic fields, the Zeeman shift 
is smaller than 200 x 0.2(1.0) kHz = 0.04(0.2) MHz. The 
crossing of the three curves yields the zero total magnetic field 
at Bappiied — 21(6)/iT (in good agreement with the expected 
earth magnetic field component of 19/iT). From this value an 
upper limit for the total residual Zeeman shift is calculated of 
0.04(3) MHz, which is negligible. 



TABLE I: Measured systematic shifts and standard uncertainty bud- 
get (la). All values in MHz. 



Effect 


Shift (MHz) 


Icr Uncertainty (MHz) 


Zeeman 


0.0 


0.0 


AC Stark repumper 


-0.4 


0.6 


AC Stark spectroscopy laser 


-0.4 


0.8 


RF Stark effect 


0.0 


1.2 


Comb calibration 


0.0 


0.2 


Statistics 




0.6 


Total 


-0.8 


1.7 



AC Stark shifts in the spectroscopy data for both the re- 
pumper laser and the probe laser were evaluated by chang- 
ing their respective intensities. As expected, both shifts were 
found to be small compared to the measurement accuracy 
(-0.4(0.6) MHz and -0.4(0.8) MHz respectively. Table ID- 
In order to assess possible trap induced shifts, the RF volt- 
age Vq was varied over a range from 50-200 V. The cor- 
responding shift was measured to be consistent with a zero 
value: A/ = 0.0(1.2) MHz. 

The '*"Ca+ 4s^Si/2 — ^P^Pi/2 transition frequency fol- 
lows from the statistical average of the measurements 
with linearly polarized light, corrected for the measured 
shifts. The total measured value adds up to / — 
755 222 765.4(0.6) + 0.4(0.6) + 0.4(0.8) + 0.0(1.2) + 
0.0(0.2) = 755 222 766.2(1.7) MHz, with the uncertain- 
ties taken in quadrature. The present result is consistent 
with the most accurate previously reported value of z^ = 
25 191.519(0.002) cm^^ (/ = 755 222 740(60) MHz). 



III. CONCLUSIONS 



We have performed an absolute frequency calibration of the 
4s ^Si/2 — 4p^Pi/2 transition in ^"Ca+ ions in a linear Paul 
trap, including a study of Zeeman and Stark shifts. The tran- 
sition was measured to be at / = 755 222 766.2(1.7) MHz, in 
a laser cooled ion cloud of T « IK. The presented value rep- 
resents more than an order of magnitude improvement over 
the previously most accurate measurement. This level of ac- 
curacy, at AA/A — 2 X 10^^, is such that for any comparison 
with state-of-the-art astrophysical data, the laboratory value 
can be considered exact. 

The employed technique for trapping and cooling ions will 
be used for spectroscopy on other ions in the near future. We 
plan to trap and sympathetically cool ions such as Mn+ and 
Ti+ with Ca+. In this way, we expect to reach improvements 
in energy level calibration of these ions, similar to the demon- 
strated accuracy for calcium. This will provide accurate lab- 
oratory wavelengths for the comparison with high redshift as- 
tronomical data in a search for a possible variation of the fine 
structure constant over cosmological time scales. 
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